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METHOD AND APPARATUS FOR SEMI-BLIND COMMUNICATION CHANNEL 

ESTIMATION 

Field of the Invention 

5 

The present invention relates to multi-carrier conununication systems, and more 
particularly to channel estimation in multi-cairier communication systems that employ 
orthogonal frequency division multiplexing (OFDM). 

10 Background of the Invention 

OFDM is a multi-carrier transmission technique, which divides available 
frequency spectrum of a communication channel into many earners, often referred to as 
sub-carriers; and adjacent sub-carriers are orthogonally phased to each other. Each of 

15 the sub-carriers is then modulated by a low rate data stream. As the sub-carriers are 
packed more closely than, for example, in frequency division multiplexing (FDMA), 
OFDM allows the frequency spectrum to be used more efficiently. In addition, OFDM 
does not require complex time switching, as in time division multiplexing (TDMA), and 
therefore does not suffer the overhead associated with time switching methods. 

20 FIG. 1 shows an OFDM system for transmitting and receiving informaUon. A 

serial stream of data symbols is provided via input 101 to a serial to parallel converter 
102 that converts the single data stream to several parallel data streams. An inverse 
discrete Fourier transform (TOFT) module 104 processes the parallel data streams and 
produces a corresponding number of orthogonal modulated sub-cairiers which arc 

25 provided to a parallel to serial converter 106. In response, the parallel to serial converter 
106 provides a serial data signal to a cyclic prefix adder 108, and the cyclic prefix adder 
108 produces a transmit data signal at output 1 10, and the transmit data signal is 
transmitted on a communication channel 

With further reference to FIG. 1, a corresponding received data signal on the 

30 communication channel is provided to a cyclic prefix remover 1 14 via input 1 12, The 
cyclic prefix remover 114 removes the cyclic prefix from the received data signal and 
outputs a single stream of data to a serial to parallel converter 116. Resultant signals 
from the outputs of the serial to parallel converter 1 16 are provided to a discrete Fourier 
transform (DFT) module 118, which provides a corresponding plurality of demodulated 

35 data streams to a parallel to serial converter 120. A serial data signal is then provided by 
the parallel to serial converter 120 to an equalizer 135 and to a channel estimator 130, 
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and the equalizer operates with the channel estimator 130 to determine the originally 
transmitted data from received data, and provide the received data via an output 135. 

The cyclic prefix is employed to address distortion in the communication channel. 
Adding the cyclic prefix comprises repeating die last few samples of each data symbol 
5 at its beginning, prior to its transmission. The length of the cyclic prefix should be 
chosen to be greater than or equal to the duration of the impulse response of the 
communication channel. This allows equalization of the channel distortion in the 
frequency domain by using a single tap scalar equalizer for each carrier, independently. 
However, in order to do this the response of the communication channel needs to be 

10 characterized. In practice only an estimate of the communication channel's 
characteristics is used, hence the need for a channel estimator. 

There are several methods of performing channel estimation, these include the 
following schemes; Pilot Symbol Assisted Modulation (PSAM), Blind Channel 
estimation, and a coded pilot method. Each of these is briefly described below. 

15 PSAM adds periodic transmissions of known symbols or pilots. Pilots comprise 

data that is known by both the transmitter and the receiver. Therefore, communicating 
pilot symbols allows the receiver to determine the difference between what was 
transmitted and what was received, and thus compensate for any variations in the 
received symbols that are caused by transmission between the transmitter and the 

20 receiver i.e. the conununicadon channel. An estimate of the characteristics of the 
communication "Channel is required to provide such compensation across time and 
frequency domains of the communication channel. When the time and frequency 
characteristics of the communication channel are varying rapidly, as in mobile 
communication applications for example, channel estimation must be performed more 

25 frequently, hence the need for more pilots to be transmitted in order to maintain reliable 
communication. Thus, reducing the available bandwidth for data transmission. 

FIG. 2 shows a PSAM scheme graphically, where both pilot symbols 202 and data 
symbols 204 are shown in a three dimensional grid across time 206 and frequency 208 
axes, and where the vertical axis represents transmission power 210 of the data and pilot 

30 symbols. In PSAM, \ho pilot symbols 202 arc inserted at intervals across time and 
frequency between the data symbols 204. Consequently, part of the signal energy and 
bandwidth of the communication channel is used for transmitting the pilot symbols 202. 
A received data signal in a PSAM scheme is passed through a 2-D Wiener filter, which 
essentially performs interpolation based on the statistics of the communication channel 

35 so as to estimate the characteristics of the communication channel between the pilot 
symbols 202, i.e. where the data symbols 204 are received. In this way, the data 
symbols 204 can be correctly recovered by taking die estimated time and frequency 
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characteristics into consideration to provide channel equalisation. One iroplemeniation 
of the PSAM scheme is in terrestrial transmission in digital video broadcasting (DYB- 
T). 

A PSAM scheme provides good channel estimation even when applied to time 
5 variant channels. However, when the normalized maximum Doppler spread is high, 
caused by fast changes in the communication channel characteristics, the frequency at 
which pilot symbols are required increases in order to track such fast changes. This 
results in more bandwidth being required for pilot symbols, up to ten percent of the 
bandwidth of the communication channel, and leaving less of the bandwidth for data 
10 traffic- 
Blind channel estimation does not use pUots. Instead, the data symbols themselves 
axe used to estimate the communication channel. Consequently, bandwidth of the 
communication channel is preserved. Several blind channel estimation schemes for 
OFDM are known, however, their tracking ability in a communication channel whose 
15 characteristics change or vary with time, Rayleigh fading time variant channels, for 
example, have not been as good as that of the PSAM scheme. 

The coded pilot method is described in US patent 5912876 by H*mimy where a 
main signal, comprising a quadrature amplitude modulated (QAM) version of a signal 
to be transmitted, and a pilot signal, are coded separately and transmitted as part of an 
20 OFDM signal. When the OFDM signal is received, the main signal portion is detected 
and an estimation of die communication channel is determined from the detected coded 
pilot signal portion. Then the detected main signal and the estimation of the 
communication channel are used to estimate the signal that was transmitted. The coded 
pilot method is simple to implement, and the coding enhances the detection of the main 
25 and pilot signals, in the consequent channel estimation process. 

However, a transceiver using the coded pilot method is necessarily more 
complicated due to the coding in the transmitter, and detection of the codes in the 
receiver. In addition, a portion of the bandwidth of the communication channel needs to 
be allocated to support the transmission of the coded signals, thus reducing the usable 
30 portion of a predetermined bandwidth. 

Hence, there is a need for a channel estimation scheme that provides good 
performance in a communication channel having varying frequency and time 
characteristics, while preserving the usable bandwidth of the communication channel. 
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Brief Suuunaiy of the Invention 

The present invention seeks to provide a method and an apparatus for semi-biind 
conimunication channel estimation, which overcomes, or at least reduces the 
abovementioned problems of the prior art. 

Accordingly, in one aspect, the present invention provides a multi-carrier 
communication system comprising: 
a transmitter comprising: 

a data input for receiving an unknown data signal for transmission; 

a known data input for receiving a known data signal; 

a signal power ratio input for receiving a power ratio signal indicating a 
ratio for combining the unknown data signal and the known data signal; and 

a data combiner coupled to the data input, the known data input and the 
signal power ratio input, the data combiner for combining the unknown data signal and 
the known data signal in accordance with the power ratio signal to produce a composite 
output signal comprising discrete data signals, wherein each discrete data signal 
comprises at least a portion of the unknown data signal and at least a portion of the 
known data signal combined in accordance with the power ratio signal, the data 
combiner having an output adapted to provide the composite output signal to a multi- 
carrier transmitter, wherein the multi-carrier transmitter transmits a transmit signal on a 
communication channel, wherein the transmit signal includes the composite output 
signal; and 

a receiver comprising: 

a multi-carrier receiver for receiving a receive signal corresponding to the 
transmit signal on the communication channel, and the multi -carrier receiver having an 
output for providing a corresponding composite signal, wherein the corresponding 
composite signal comprises corresponding discrete data signals, and the corresponding 
composite signal being shaped by at least one signal shaping characteristic of the 
communication channel; 

a channel estimator having a known data input for receiving the known data 
signal, an input coupled to receive the corresponding composite signal, an input coupled 
to receive the power ratio signal, and an input for receiving at least one estimate of the 
unknown data signal, the channel estimator for estimating the at least one signal shaping 
characteristic of the communication channel from at least the corresponding composite 
signal, the at least the portion of the known data signal of at least some of the 
corresponding discrete data signals, the power ratio signal and the at least one estimate 



of the unknowTi data signal, and the channel estiniator having an output for providing at 
least one estimated communicaiion channel characterising signal; and 

an equalizer coupled to receive the corresponding composite signal, the 
known data signal, the power ratio signal and the at least one estimated communication 
channel characterizing signal, the equalizer for configuring at least one of its signal 
shaping characteristics to compensate for the at least one signal shaping characteristic of 
the communication channel, the configured equalizer for shaping the corresponding 
composite signal accordingly, and the equalizer having an output for providing at least 
one subsequent estimate of the unknown data signal. , 

In another aspect the present invention provides a multi-carrier transmitting 
system comprising: 

a data input for receiving an unknown data signal for transmission; 

a known data input for receiving a known data signal; 

a signal power ratio input for receiving a power ratio signal indicating a ratio for 
combining the unknown data signal and the known data signal; and 

a data combiner coupled to the data input, the known data input and the signal 
power ratio input, the data combiner for combining the unknown data signal and the 
known data signal in accordance with the power ratio signal to produce a composite 
output signal comprising discrete data signals, wherein each discrete data signal 
comprises at least a portion of the unknown data signal and at least a portion of the 
known data signal combined in accordance with the power ratio signal, the data 
combiner having an output adapted to provide the composite output signal to a multi- 
carrier transmitter, wherein the multi-carrier transmitter transmits a transmit signal on a 
conmiunication channel, wherein the ttansmit signal includes the composite output 
signal. 

In yet another aspect the present invention provides a multi-carrier receiving system 
comprising: 

a multi -carrier receiver for receiving a receive signal on a communication 
channel, where in the receive signal includes a composite signal, and the multi-carrier 
receiver having an output for providing the composite signal, wherein the composite 
signal comprises discrete data signals, and wherein each discrete data signal comprises 
at least a portion of an unknown data signal and at least a portion of a known data signal 
combined in accordance with a signal power ratio signal, the composite signal being 
shaped by at least one signal shaping characteristic of the communication channel; 

a channel estimator having a known data input for receiving the known data 
signal, an input coupled to receive the composite signal, an input coupled to receive the 
power ratio signal, and an input for receiving at least one estimate of the unknown data 



signal, the channel estimator for estimating the at least one signal shaping characteristic 
of the communication channel from at least the composite signal, the at least the portion 
of the known data signal of at least some of the discrete data signals, the power ratio 
signal and the at least one estimate of the unknown data signal, and the channel 
estimator having an output for providing at least one estimated communication channel 
characterising signal; and 

an equalizer coupled to receive the composite signal, the known data signal, the 
power ratio signal and the at least one estimated communication channel characterizing 
signal, the equalizer for configuring at least one of its signal shaping characteristics to 
compensate for the at least one signal shaping characteristic of the communication 
channel, the configured equalizer for shaping the composite signal accordingly, and the 
equalizer having an output for providing at least one subsequent estimate of the 
unknown data signal. 

In still another aspect the present invention provides a method for determining 
received data in a multi-carrier communication system, wherein a received signal 
includes a composite signal received on a communication channel having transmission 
characteristics, wherein the composite signal comprises a plurality of discrete data 
signals spaced in lime and frequency, and wherein each discrete data signal comprises a 
data portion and a pilot portion, wherein the data portion comprises one of a 
predetermined group of symbols, the method comprising the steps of: 

a) defining ia set of the plurality of discrete data signals; 

b) setting a predetermined number of iterations; 

c) receiving the composite signal on the communication channel; 

d) selecting one of the plurality of discrete data signals to be estimated; 

e) selecting a group of the plurality of discrete data signals, wherein each discrete data 
signal of the group of the plurality of discrete data signals is relevant to determining the 
one of the plurality of data signals; 

f) normalising the group of the plurality of discrete data signals using at least the pilot 
portion of the one of the plurality of discrete data signals; 

g) estimating the transmission characteristics of the communication channel using the 
normalised group of the plurality of discrete data signals in step (f); and 

h) estimating the data portion of the one of the plurality of discrete data signals using 
the communication channel having the transmission characteristics estimated in step(g) 
and at least the pilot portion of the one of the plurality of discrete data signals. 
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Brief Description of the Drawings 

An embodiment of the present invention will now be more fully described, by 
way of example, with reference to the drawings of which: 
FIG, 1 shows a prior art communication system; 

FIG. 2 shows a graphical representation of data and pilot symbols communicated 
on a communication channel of the communication system in FIG. 1; 

FIG. 3 shows a communication system in accordance with the present invention; 

FIG. 4 shows a graphical representation of data and pilot symbols communication 
on a communication channel of the conmiunication system in FIG. 3; 

FIG. 5 shows a flow chart diagram of a process in a receiver of the 
communication system in FIG. 3; 

FIG. 6 shows a graphical representation of the APSB scheme in the 
communication system in FIG. 3; and 

FIG, 7 and 8 shows comparative simulation results of the performance of the 
communication system in FIG.3 and prior art communication systems. 

Detail Description of the Drawings 

The present invention combines the advantages of the PS AM and the blind 
schemes to produce a scheme where data and pilot symbols arc combined prior to 
transmission, and separated when received at a receiver. In addition, separation of the 
pilot symbols at the receiver is accomplished by treating the data as noise and applying 
an iterative process to detect die data symbols. An embodiment of the present invention 
will now be described. 

In FIG. 3 a conmiunication system 300 that communicates data using an added 
pilot semi-bhnd (APSB) scheme, in accordance with the present invention, comprises a 
combiner 305 having a data input 307 for receiving a data signal 306, a pilot input 309 
for receiving a pilot signal 308, and an input 311 for receiving a data-to-pilot power 
ratio (DPR) signal. The data signal 306 comprises a series of data symbols while the 
pilot signal comprises a series of pilot symbols. The DPR signal includes information 
on power levels of the data signal 306 and the pilot signal 308 that are to be combined. 

The combiner 305 includes a data conditioner 312 that receives the data signal 
306 and the DPR signal and provides a conditioned data signal having a power level as 
indicated by the DPR signal. Similarly, the combiner 305 includes a pilot conditioner 
313 that receives the pilot signal 308 and the DPR signal, and provides a conditioned 
pilot signal having a power level as indicated by the DPR signal. 
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The combiner 305 also includes an adder 314 for combiniiig the conditioned data 
signal from the data signal conditioner 312 and the conditioned pilot signal from the 
pilot signal conditioner 313, by combining a series of conditioned data symbols and a 
series of conditioned pilot symbols, and providing data-pilot signals 310 via an output 
5 318 of the combiner 305. The data-pilot signals 310 comprises a series of discrete data 
signals, v^/here each discrete data signal has a data portion and a pilot data portion, and 
where the power levels of the data portion and the pilot portion arc in accordance with 
the DPR signal. The data portion represents at least one data symbol 

The data-pilot signals 310 are then processed, in sequence, by the serial to parallel 
10 converter 102, the EFFT 104, the parallel to serial converter 106, the cyclic prefix adder 
108, as described earlier, and a transmit data-pilot signal from output 310 is transmitted 
on a communication channel. 

A corresponding receive data-pilot signal 314 received from the communication 
channel at input 315, is processed sequentially by the cyclic prefix remover 1 14, the 

15 serial to parallel converter 116, the discrete Fourier transform 118 and the parallel to 

serial converter 120, as described earlier. An output data-pilot signal 316 from output - 
317 of the parallel to serial converter 120 is then provided to an APSB equalizer 325 
and to an APSB channel estimator 320. The APSB chaimel estimator 320 and the APSB 
channel equaliser 325» each include an input 319 for receiving the pilot signal 308, and 

20 each include an input 3 1 8 for receiving the DPR signal. The pilot signal 308 and the 
DPR signal rriay be stored in a memory (not shown). The APSB equaliser 325 operates 
with the APSB channel estimator 320 to determine the originally transmitted data 306 
from the data-pilot signal 316 using the pilot signal 308 and the DPR signal, and 
provides the originally transmitted data 306 via output 330 of the APSB equaliser 325. 

25 This is accomplished with an iterative process where a channel estimate h(n) is made 
by the APSB channel estimator 320, and a data estimate is then made using the 
previously obtained channel estimate and provided via output 330. The output 330 is 
coupled to the APSB channel estimator 320 to provide the channel estimate for a 
subsequent iteration of channel estimation and data estimation. When a predetermined 

30 number of iterations have been perfonned, the originally transmitted data 306 is 
determined and provided via output 330 of the APSB equaliser 325. 

With additional reference to FIG. 4, which shows the APSB scheme graphically, 
symbols of the data-pilot signal 317 are shown as discrete signals 405 in the three 
dimensional grid across time 206 and frequency 208 axes, and where the vertical axis 

35 represents transmission power 210. Each discrete signal 405 comprises two portions, a 
data portion 407 and a pilot portion 409. The proportions are defined by the DPR signal, 
which is. as defined earlier, the ratio of the data signal power to the pilot signal power. 
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There is an optimum DPR where the performance of the communication system 300 
would be best. 

In accordance with the present invention as described, the pilot is transmitted 
along with the data simultaneously and in the same frequency band, therefore 
5 transmission of the pilot advantageously does not consume bandwidth of the 
communication channel. 

The data-pilot signal 314 received from the conmiunication channel contains 
information about the communication channel, and the pilot portion of the data-pilot 
signal is known. Then, by treating the data portion as noise, an initial estimation of the 

10 communication channel can be determined, and a data estimate obtained, with the initial 
estimation of the communication channel, using a least squares approach. With partial 
knowledge of the characteristics of the communication channel via the channel 
estimates and the data estimates, further iterations of channel and data estimation are 
performed, thereby improving the accuracy of the estimation and allowing an accurate 

15 estimate of the transmitted data to be made after a predetennined number of iterations. 

With additional reference to FIG. 5, a process 500 for estimating the transmitted 
data starts 505 with selecting 510 a set of cells from the received set of data-pilot 
cells S. Next the predetermined number of iterations is set 515, and a first cell, for 
example cell 430 at (n,l) is selected 535. Subsequently, cells relevant to the estimation 

20 of the cell 430 are selected 540. Here, the selected cells are the eight cells adjacent to 
the cell 430 as typically the closest surrounding cells are selected. Using the time 206 n 
and frequency 208 / axes as references, the cell 430 is identified as cell (n,l) and the 
relevant selected cells as S^(nJ), All the relevant selected cells S^(nJ) are then 
normalised 542. for example by dividing the selected relevant cells by the pilot symbol 

25 portion 409 to allow a suitable channel estimate to be made. An initial channel estimate 
545 is then made using a 2-D Wiener filter, where the data portions 407 of the selected 
relevant cells are treated as noise, and the result is an estimate of the communication 
channel at (nj). Subsequently, a data estimate 550 is made using the least squares 
approach, and a decision 555 is made as to which symbol was transmitted. The decision 

30 Step 555 may involve processes such as interleaving and coding before a decision is 
made. 

A determination 560 is then made as to whether all the cells in the set S have hctn 
determined. When not all the cells in the set 5 have been deienmned, then the next cell 
of the set S of cells is selected 570, and the process 500 returns to step 540 of selecting 
35 relevant cells to the newly selected cell. The looping back through step 570 continues 
until all the cells in the set S have been determined. When ajl the cells in the set S have 
been determined, the determination at step 560 is true, a counter (not shown) indicating 
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the number of iterations i is incremented 562, and a deteimination 565 is made as to 
whether the number of iteration i have reached the predetermined number of iterations 
set in step 515. When the number of iterations i has reached the predetermined number 
of iterations set in step 515, the process 500 ends 567. However, when the number of 
iterations / has not reached the predetermined number, the process 500 returns to step 
535 and repeats as described above. 

Returning now to FIG. 3 a functional description of ^the communication system 
300 now follows. The rules for notation that are employed are: 

in the time domain, signals at index (n,/) are written* as h(nj), vectors as h(n), and 
matrices as H(n); 

in the frequency domain a subscript F is added (e.g. hffn,/j, tuin), Hffnjj; ^nd 
sub- vectors/matrices are capped with a tilde e.g. h or H , and estimated entities 
with a hat e.g. h; 

superscripts ^ ^ and * denote the operations matrix transpose, Hermitian and 

conjugation, respectively, and E denotes expectation; and 

elements of vectors/matrices are denoted as [h]/[H]-. 
. The elementary sample period of an OFDM system model is T, the number of 
subcarriers is L, the number of cyclic prefix samples is Z?, and the maximum number of 
channel response samples is 



and. the OFDM symbol period is defined as 

The complex baseband representation of the communication channel 605, 
particularly a mobile wireless communication channel, impulse response at time t is 
described by 



5+1 (<D) 

while the total number of samples in one OFDM symbol is 




(1) 



where T/fr) and Y/(0 arc the delay and complex ampJitude;of the /"* path, respectively. 
The power delay profile of the channel is defined as 



p(X)^EmT)h\t,T)] ; : (2) 



Assuming an exponentially decaying power delayjprofile with 

J 

p(T)=Acxp(-T/T^) I 
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where is a parameter of the channel known as the too 



mean squared delay spread, 



and Aish normalizing constant. The normalized T„^ is given as x^r. 

Due to the relative mobility of the transmitter and receiver, a maximum Doppler 
spread. /o, will occur in the received signal. This is accounted for with a time-variant 
channel where a high^ implies a fast varying channel. The normalized maximum 
Doppler spread is defined as/^r,, based on the assumption that the channel coefficients 
are time invariant over each OFDM symbol period T,. 

A function which is useful for analysis in OFDM ils the time-variant transfer 
function obtained from the Fourier transform of equation 
delay T, which produces the equation below. 



;i) above, with respect to 



h{tX) exp(-/2Jl/c) d% 



(3) 



Assuming a wide sense stationary uncorrelated scattering (WSSUS), and a 
15 Rayleigh fading channel with Jakes* spectrum, the autocorrelation of the channel is 



separable in time r, and frequency/, and can be written as 

where At = t- f -ind Af= //'.We have 

l-exp(-Z)r(l/T,^ +;2/!A/)) 



follows. 



C4A) 



(1 - exp(-Dr/T_))(l + jlmj) 
1 



l + jlAfT, 



(4B) 



assuming [(D)/(W7}] » I, and 

r,(Ar)=yo(27ifbA/) (4C) 

with /(,(■) being the zcroth order Besscl function of the first kind. 
The discrete channel model of a communication channel ijf order B( < D) can be 
described as /i(«,0 = h(t = nT,,T = IT), assuming that Ti^(t\)is uniformly spaced at 
intervals of T and that P(^^)=0 for t ) t,^ = 57 and represented as a vector at time 
nas h(rt)=[/»(n,0)...A(n,L-l)f 

=[h(n.0)...h{n.B)0...0f eq^^^ 
where C is a complex number. 



(5) 
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The channel coefficients in the frequency domain are obtained in a similar way 
to equation (2). but in discrete time as . 

hp(n)=VZFh {n)=[hp{n,Oy..hpin,L-l)Y 



(6) 



5 where F is the L X L unitary discrete Fourier transfonn (DjFT) matrix with 



^^^P(~^^'^ £-)/ VI, and F^is the corresponding inverscDFT (IDFT) matri?c. 
The auto-correlation of the discrete channel, \fi^(7i,l), can be obtained from the 
equations (4A), (48) and (4C) by replacing ^f=I/(LT) apd At = nTs, 

The received OFDM signal may be considered a digital signaJ in a 2-D cell 
structure with indices (n,l)ES- 

where S={{nj):0 <n<N -l,0<l<L-i} , and where is die total number of cells in 
the lime direction, n, and L is the number of sub-carriers as defined earlier. 

A model for a 2-D Wiener filter will now be described, as such a filter is used to 
form estimates from sampled signals at the receiver. Estimates will be obtained at 
index inJ)sSp from the sampled signals at indices (n\F\)eSi^'^\ where 5^ and 5^'*''> 
are subsets of 5, and where Sp =u„ ^ S^p'^K The definitions of Sp.Sp and S^p'^^ be 
provided later. i 

In a 2-D wide sense stationary (WSS) stochastic ptocess hf(nj) 
which contains information of the desired signal hp>(nj)^^d is corrupted by coirelated 
noise hp(nj)vp(nj) and additive white Gaussian noise (AWGN), Up(nJ), as follows: 



hpinj} - hp(nJ)-\-hp(nJ)vp{nJ)-\-Up(nJ) 



(7) 
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The notation hp(nj) and indicate that the channel frequency response 

is estimated from its noisy samples. It is assumed that Vp{n, I) and Uf:{n,l) are 
white, and that /z/:'(n,/)»v^(w,/)andu^(n,/)are mutually u;i-correlated, and zero mean 
stochastic processes. In the prior art PS AM scheme, the correlated noise component is 
absent, in contrast, for the APSB scheme of the present injvention, as described, the 
correlated noise component is advantageously reduced through successive iterative 
measures. 

For the 2-D Wiener filter, the estimator for A^(n,Z)is given as: 



(8) 



where w(n,f;nT ) is the weight applied on hp(ri J){o estimate hp{nJl 

Based on the assumption that a total of Ay^/ and A/ pilots are used in the n and / 
direction, respectively, for each estimation of hp{nj) , a \ector fi(n,/)eC^^^^ ^ ^ is 
performed by stacking the elements from hfin' ,/* ), V(n* , T )sS^'^] and a vector 
w (n,/)eC^^^^^ Is formed by stacking the conjugate elements of >v(n,/;n'r). 



13 



Equation (8) can be re-written as follows. 

To minimize the difference hp{nj) — h(nyl) in the mcani i 
optimum tap-weight vector, the orthogonality principle is now applied to equation (9), 
5 resulting in the equation below. 
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35 



(9) 

square sense and obtain the 



E [lihp(n,l)-hp(n,l)).h}(n\n] = 0, v{n".r} e 4"''^ (lO) 
Substituting equation (8) into equation (10), the Wiener-Hopf equation, with optimum 
tap weights Wo(n,i;n' ,f ) , is obtained as provided below 

"2 Wo(nJ;ri J)E[hpin' ,r)h;(n\n] = Mh.Mh'^in-'J")] (H) 

Defining the cross-correlation of hp{n,l)andhpfyt,l) and the autocorrelation 
of hp(n,l), respectively as 



.n] 



(12) 



(13) 



Letting. r-^^^(n,0=£[ h pin,l)h'ir{n,l)]eC^''^'''^ be the cross- 
correlation vector and R^^^ (n,0=£[ h F{n,l) hpinJ^jeC^^'^'-''^''^'- be the 
autocorrelation matrix hp(n,l) formed from the elements! of equations (12) and (13). 
respectively. Therefore, the Wiener-Hopf equation can be re-writtcn in matrix notation 
25 as 

Wo(«,/)=Ri-^r^^./n.O (14) 

Assuming that hf(nj),v{n,l) and u{n,l) are mutually un-correlated and wide- 
sense stationary white stochastic processes, using the definition provided by equation 
30 (7), and letting An and AZ be the discrete time and frequency difference indices, 
respectively, equations (12) and (13) can be written as follows. 



%hf = (^ A/) + r^^hr ^0 X 

+ r^^v^(An,A/) 

/),^,,^(An,AO 



= { 



(15) 



forAn,AZ = 0 
otherwise 



(16) 



i 

i 
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where r/,^^^(An,A/) is the autocorrelation matrix ft^. 

R ^- ^ and ^ (n,/) can be simplified as * 

And, where R hphp is the NL by autocorrelation matrix of hp, then R/.^;,^can be 
obtained from equation (4). If R andal are known, the optimum weight is 



(20) 



(19) 

10 The mean square error (MSE) can then be obtained as 



4K(«»0-^^(n,Z)p]v(n,065o^ 
The minimum MSE (MMSE) at index (n,l) is then 

equation (9) into equation (20). and using the optimum weight factor in equation (19)» 

AfM5£(n,/)=(T^-w^(n,/)r;,^;^/n,0 (21) 



obtained by substituting 
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For simplicity, or/j=l for all subsequent equations. 
Now letting the transmitted data in the frequency |omain be 
df{n,l)y{nJ)eSQ , where n is the discrete time index |and / the discrete frequency 
index. In an OFDM system, n refers to the OFDM block index, while / refers to the sub- 
20 earner index, and pilots are denoted as /7;r(n,Z), V(n./)e sl>. It is assumed that Pf(nJ) is 
deteraninistic and selected from a fixed set of alphabets, \yhile dj^inj), the data is a 
zero mean stochastic process. The variance of ppirt.l) is ^denoted as s^o^ , and in 
practical implementations £^((1 . | 

In accordance with the present invention, as described herein, pilot signals and 

r 

23 data signals arc added together and co-cxist at all time, n,]and frequency points , /, in 
order to conserve bandwidth, that is, S^^Sp =5. The signal after the data signal and the 
pilot signal are combined is defined below. , 

Xf{n,r)^dp{n,l)+p^{nJ\ ^M^S (22) 

and cr^ = (l+£:^)o-3. Other definitions follow. I 
30 Data -to- Pilot Power Ratio iDPR)=al /(ej o'2)=l/ 

Signal -to- Pilot Power Ratio(SPR)=a^ /(i^crj)= I-f 1/ £^ = H-DP/? 
Assuming that the orthogonality of the OFDM system 600 is maintained, the 
signal after OFDM demodulation is 

35 =hf(nMdf(nJ)-hpjrinJ))ApinJ) (23) 

where ft^(n,Z)is an additive Gaussian noise (AWGN) stochastic process. 



# 
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In accordance with the present invention, as described herein:, channel-data 
estimation from the received signal yf(n,/) is performed injan iterative manner. For ease 
of description, the first iteration and subsequent iterations jWill be described separately. 

FIRST ITERATION | 

Channel estimator 320 is a 2-D Wiener filter, as characterized earlier, estimates the 
channel response. First, the received signal yf(n,l) in equation (23) is normalized in the 
first iteration, thus. 

=/2;.(rt,0+^F(nJ)i/r(".0/Pf(nJ)+Z7^(n.O/p;r(n,/) 
=/z;,(n,/)+/2p(n./)v^i(7j,/) + u[, (7i,Z) (24) 



where Vp^(nJ) = df{nj)/ Pf{n,t) is a data-dependent noise introduced due to the addition 
of the pilots to die data, and u^^{nj) = b^{nj) I ppin.t) is A^GN. The subscript "1" 
indicates that the notation is specific to the first iteration. Similar use of subscript "r" 
will be employed for the iteration. 

Assuming that the pilots, Pfin^l), are selected from] ; 
symbols, we obtain 



and similarly 



a set of constant modulus 



■ ^5 _ 



=^DPR 



2 _ Oh _^^^s?R 



sr^R 



(25) 



(26) 



where, as defined earlier. SPR=(l + e^)/e^,andSNR=cr^ 



the optimum tap weights can 



Equation (24) is in a similar fonn as equation (7) and 1 
be obtained in a similar form as in equation (14) using time A^^ and frequency A^, 
sampled signals. In this case, the autocorrelation matrix R^^^ » is independent of (n.l), 
however the cross correlation vector r^^;^^ (".0 is usually dependent of (n,l) . 

FIG. 6 shows a graphical representation in the time/frequency domain, and ^p'^^ 
is selected such that (n,l) is in the center of S^p'^\ The set o|f indices encompased by 
will change, and may be visualized as a sliding window, as fn,/).changes. The relative 
index of (nj) and 5p"''* however, will not change. Assuming wide sense stationarity, 
R i^^f^^i^yl) will also not change. Next, for example, the cyoss-correlation vector 
''vv^"»'^ when estimating index 601 and 602 would be the same. When (n,/) is near to 
the boundary of 5/>, such as when estimating index 603, the relatave indices of (n.O and 
V*'* would be different. Thus,r/;^^^(«,/) would be differejnt in this case. 
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Now, the MMSE estimator for the first iteration fojllows 

^Fi(''.o=w s;,(/i./)hn(«./) 

Thus, the first esumate of JC^n,/) is 
Then, in view of equation (22), the estimate for dp{n,l) is obtained as 



(27) 



(28) 



_yF{n,l). 



_hp(^n,lXdFin,l)+pp{n,l))+bp{n,l) 



10 



hn(n,l) 



If the channel estimation is perfect, i.e. ^^(n,?)= hp{n,l). 



Then, a decision device or slicer, as is known in the art, is 
an estimate dfr^{n,l) of dp(n,l). 



15 



-p,(.nj) (29) 



then equation (26) becomes 



I 



(30) 



used on dir^(n,l) to obtain 



SUBSEQUENT TTERATTONS 
For the second iteration, the normalization is carried out using 5„(n,/) + pp(n,l). 
The normalization tries to remove the data-dependent 
the first iteration as characterized in equation (24). The normalization for the second 
20 iteration is carried out as follows: 



hF-,(n,l)=^ 



)'f("./) 



dpi(n,0+pir(n.l) 
_ hr(n, l)(dp(n, I) + p^(n, /)) } b^(n,l) 



25 



dp^in,l)+pj,in,l) 



(31) 



Making the substitution ^p,(n,l) = dp{n,l) which is a good approximation when 
the probability of symbol error is small, equation (31 )becomes 



dpiinJHp 



30 where 



«rt(".0= 



(32) 



(dnin,l)+prin.l)) 
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Equation (32) has the form of equation (7) with crjj ~^ ^P^^ v;;2(n,/)=0. 

Then based on the assumption that E^{{l,df(nJ) is modulated by quadrature phase 
shift keying (QPSK), and PF{nJ)e[eoa(±l±j)/'sf2} , a QPSK consteUation 
\Sfiin,l)+pp{nJ)^c^ have four possible values, each being equally likely to occur: 

€r^\(li-j)'^e(±l±jf/2 . It should be noted that only tie first quadrant of the QPSK 
signal constellation has been considered due to its 2-D symmetry. Hence, the required 
variance reduces to = ^2 E\\d,,M+p,in.i f] 



=(!+£*) 



(l-2e^+e<)(jl+£V^ 



=(l+£'')(l+2£^+-)SM?"' 
This approximation can be made when £^ ({1. 



(33) 



equation (9) and (14) by 



The 2-D Wiener filtering is applied accordingly in 
setting by setting ^2 =^0-^2= 0 ^^'^^ "° correlated noise component) and 

~^ui Then the channel estimator is 



Finally, the estimate ^^jC^^Ois obtained by following the steps as in equations 
(28) and (29). 

The subsequent iterations can be extended from e second iteration. However, 
the assumption dp^in.Vj-dpinJ) made in equation (32) viU be more accurate in 
25 subsequent iterations and will therefore result in a better estimate for 



(34) 



hf^{n,l),hp^{nj), etc. As more iterations are perFormed, 
converges to the actual channel response 

In the prior art PSAM scheme, for every one pilot 



(^yv^L data cells. The effective average signal po vcr after pilot insertion is 



scheme. 



y?^ + (<?A-O l (35) 
where rj^ > 1 is the ratio of the power of the pilot lo thai of the data for PSAM 



Lhe channel estimate 



cell, there are on average 
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The signal-to-noise ratio (SNR) of the PSAM scheme can 
expressed as 



# 



thus be 



SNR^al lal is the SNR of the system without the presence of pilots, while SNR; 



5 represents the additional SNR required to compensate for 



been allocated to the pilot. Therefore, it follows from equations (35) and (36) that 

SNR,,^,=mog,^{alla])dE 
=101ogio 



Since part of the bandwidth is used for pilot transijiissiOD, there is a reduction in 
the effective data rate. The percentage bandwidth loss is 



on (36) for proper comparison 
additional power is incurred 



When power boosted pilots {i:e. r;)!) are used, we 
actual SNR would have to be adjusted according to equati 
with other schemes. In the case when t/= 1 , however, no 
under the definition provided by equation (37). Regardlesk of the value of rj, a 
bandwidth loss is still incurred as seen from equation (38 

For the APSB scheme, in accordance with the preient invention, as described, 
the SNR is defined in a similar way as in equation (36). h >wever with the APSB 
scheme 



(36) 



the loss of power which has 



(37) 



(38) 



note that SNR,^^ > 0 dB. The 



( SPR 
\SPR-l 



^101og,nl^:^^|rf5 



(39) 



In contrast to the prior art PSAM scheme, the APSB scheme SNR,^ is always 
greater than 0 dB since £)0. On the other hand, although the prior art PSAM scheme 
suffers from bandwidth loss, the APSB scheme advantageously has Wj^^ = 0. 

Both the prior art PSAM and the APSB schemes require sufficient OFDM 
symbols to be received before an optimum channel estimate can be obtained. 
Depending on the selection of the indices of the sampled signal, Sp"-^, different 
processing delays would be incurred in order for an optinium estimate to be made. 



1 
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When a non-causal 2-b Wiener filter is used, such that the! index (n,/) is at the center of 
Sp^"-^, the processing delay is then Sfj{AN- 1)/2 for the prior art PSAM scheme. For 
Tjie APSB scheme, the processing delay is the same with = 1 . 

The results of simulations will now be presented. For the simulations, the 
common system parameters were set as L = 32, D =8, B =7, '^rms IT- 0.5, 'J'm^ /T = 8 
and SNR = 20dB. It was assumed that the SNR is known ^nd that (J,^= 1, and a QPSK 

signal was used. In addition, 
of the 2-D indices of Sp^""'^, 
estimate the channel at 



/2or/>L-(A^+l)/2 

out. For the other indices, an 



signal constellation for the transmitted 'signal and the pilot ; 
5p ""^ was selected such that the index (n,0 is in the center 
This results in the best MSE performance. It is possible to i 
10 indices 

rt<(A^ -l)/2orn>A^-(Ay^+l)/2, and /<(A£^-1)1/ 
in such a mariner where a ccntered-inteipolation is carried i 
off-centered interpolation is required to estimate the channel, giving some degradation 
of the MSE. 

15 With reference to FIG. 7, consider first, the bit error rate (BER) performance at a 

maximum normalized Dopplcr spread,;/or,. of 0.005. Sett ng A/^]= A^ = 15, the BER for 
the APSB scheme using different e(SPR) for the first to ti e fifth iterations. A high SPR 
means that low amount of power is used for pilot transmission while a small SPR means 
that low amount of power is used for data. Both extreme cases imply that hf(nj) will 

20 be badly estimated. Hence, for each iteration, there exist an optimum point whereby the 
SPR will mininiize the BER. It is obseSrved that as the nuriber of iterations increases, 
the performance improves, but at a diminishing marginal rniount. 

For the prior art PSAM scheme, as we increase TJ, the channel estimation 

improves. However, since SNRio„ also] increases, less povs er is allocated for the data for 
25 a given SNR. Thus, the improvement of the BER brought about by better channel 
estimation would be offset at some point as T) increases. Gonsequently, there is an 

optimum T] that minimizes the BER as well. 

I 

With both the prior art PSAM scheme and the APSB scheme, performance 
improves when larger filter taps are us^ed i.e. larger A;^ and A^. For the PSAM scheme, 

30 the performance also improves for smaller A^^ and A^ when closer sampling intervals 
are employed, while for the APSB scl^eme, more iterations (K) results in better 
performance, although the marginal performance gain would decrease for both cases. 
Hence, the system design is flexible since increased complexity of the PSAM scheme, 
and the APSB scheme would improve|performancc. However, the APSB scheme, in 

35 accordance with the present inventionias described, advantageously does not suffer 



I 

f 
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bandwidth loss in contrast to the PS AM scheme which wil 
bandwidth. 



In general the BER is relativelyirobust to the selection of E and rj, and is 



consistent for other values of A^^ ( = AJ, A// ( = A J, and di 
spread. Optimum values selected are as 6^ = 0.43 and Tl^, = 

With reference to PIG. 8, assuming that the SNR is 
plot of BER vs SNR for the prior art switched based blind 



always, have loss of 



fcrent values of Doppler 
4/3 that minimize BER. 

known, the graph shows a 
scheme 801, andPSAM 802 



scheme, and the APSB scheme 803 and 804, using 6^ - 0.43 and = 4/3, at^r, = 
0.005. It is observed that when A^y = A^.= 1 5» plot 803. the SNR of the of the APSB 
10 scheme is 1 dB worse than the PSAM scheme. However, the BER performance of the 
APSB scheme is superior to the switched based Wind scheme. 

The present invention, as described, provides an added pilot semi-blind scheme 
that does not consume bandwidth, and is suitable for use in a mobile communication 
system. ; 
IS This is accomplished by adding data and pilot data at the transmitter prior to 

transmission on a communication channel, and using a 2-ID Wiener filter to recover the 
transmission characteristics of the communication channe using the pilot data. Then 
through an iterative process of estimating the communication channel and estimating 
the data, the communication channel can be estimated and the transmitted data 
20 recovered. 

The present invention provides a method and an apparatus for semi-blind 
conmiunication channel estimation, which overcomes, or ; 
abovementioned problems of the prior art 

It will be appreciated that although only one particlilar embodiment of the 
25 invention has been described in detail, various modificaticns and improvements can be 
made by a person skilled in the art without departing from 1 
invention. 

We claim: 



iit least reduces the 



the scope of the present 



